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It is shown that a Coulomb suppression of the stellar enhancement factor occurs in many endother- 
mic reactions at and far from stability. Contrary to common assumptions, reaction measurements 
for astrophysics with minimal impact of stellar enhancement should be preferably performed for 
those reactions instead of their reverses, despite of their negative Q value. As a demonstration, the 
cross section of the astrophysically relevant *^Rb(p,n)**^Sr reaction has been measured by activation 
between 2.16 < iJc.m. < 3.96 MeV and the astrophysical reaction rates at p process temperatures 
for (p,n) as well as (n,p) are directly inferred from the data. Additionally, our results confirm a 
previously derived modification of a global optical proton potential. The presented arguments are 
also relevant for other a- and proton-induced reactions in the p, rp, and vp processes. 

PACS numbers: 26.50.-l-x Nuclear physics aspects of novae, supernovae, and other explosive environments, 
24.60.Dr Statistical compound-nucleus reactions, 27.50.-|-e 59 < A < 89 
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I. INTRODUCTION 

Astrophysical reaction rates are central to tracing 
changes in the abundances of nuclei by nuclear reactions. 
They provide the temperature- and density-dependent 
coefficients entering reaction networks, the large sets of 
coupled differential equations required to study nucle- 
osynthesis and energy generation in astrophysical envi- 
ronments. The reaction rates are computed from reac- 
tion cross sections which, in turn, may be predicted in 
theoretical models or extracted from experiments. In ad- 
dition to the difficulties arising in the determination of 
the cross sections, the conversion to reaction rates is fur- 
ther complicated by modifications of the rates in a hot 
plasma and the fact that the rates of forward and reverse 
rate for the same reaction have to be consistent to ensure 
numerical stability and proper equilibrium abundances at 
high temperature. Both issues can be addressed at once 
by accounting for the thermal population of target states 
which leads to stellar rates obeying a reciprocity relation 
between forward and reverse rate. Using this reciprocity, 
knowledge of the rate in only one direction is needed 
because the other reaction direction can be directly com- 
puted from that rate, thus ensuring consistency. 

For numerical reasons, further elaborated in Sec. IIIB( 
it is usually preferable to start from the rate of a reaction 
with positive reaction Q value when computing the rate 
for its inverse reaction. Even more importantly, experi- 
mentalists want to determine rates as close as possible to 
the actual stellar rates, i.e. rates with minimal thermal 
population effects of the target. Again, it can be ar- 
gued that this is the case for exothermic reactions. This 
led to the commonly applied rule that measurements of 
exothermic reactions are more important than those of 
endothermic ones. In this paper we show that there is a 



considerable number of reactions for which a suppression 
effect brings the stellar rate of an endothermic reaction 
closer to the laboratory value than its exothermic coun- 
terpart. 

As an example of how to exploit this suppression ef- 
fect and to obtain stellar rates from a measurement of 
an endothermic reaction, we experimentally studied the 
reaction *5Rb(p,n)^5Sr, having Q = -1.847 MeV. The 
importance of the reaction is manifold. In the last sev- 
eral years a number of proton capture cross section mea- 
surements with relevance for 7 process studies have been 
carried out (see, e.g., [l| and references therein). The 
7 process was shown to synthesize p nuclides (proton- 
rich isotopes not accessible to the s and r processes) by 
a series of photodisintegrations of stable nuclides in hot 
layers of massive stars [I, [3, 0, [1| • Recently, systematic 7 
process simulations found not only that photodisintegra- 
tion reactions are important but also that (p,n) reactions, 
and in particular *^Rb(p,n)*^Sr, strongly influence the fi- 
nal p abundances @ . Additionally, this reaction is well 
suited to test the optical potential used for calculating 
the interaction between protons and target nuclei. 

We commence by outlining the theoretical background 
regarding stellar rates and the suppression effect in Sees. 
Ill Al and lllBI The results of a large-scale study of the ef- 
fect in the full extension of the nuclear chart are discussed 
in Sec. Ill CI Focusing on the reaction *^Rb(p,n)*^Sr, 
their relevance is discussed in Sec. lIII K[ the experimental 
details are provided in Sees. IIII Bl HITdI and the astro- 
physically relevant rates are derived in Sec. IIII Gl Ad- 
ditionally, Sec. IIII Fl discusses implications of our new 
experimental results for the proton optical potential. Fi- 
nally, a summary is given in Sec. IIVI 
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A brief account of our findings was already given in 
Q- The present follow-up paper expands the discussion 
and also provides additional results in all parts of this 
investigation. 
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FIG. 1: Schematic view of the transitions (full arrows de- 
note particle transitions, dashed arrows are 7 transitions) in 
a compound reaction involving the nuclei A and F, and pro- 
ceeding via a compound state (horizontal dashed line) with 
spin Jc and parity ttc in the compound nucleus C. The reac- 
tion Q values for the capture reaction (Qcap) and the reaction 
A— >F (Qp) are given by the mass differences of the involved 
nuclei. The effective cross section g^^ (Eqs. |4] and [6]) for a 
reaction type is a sum over all energetically possible transi- 
tions to bound states (capture: in nuclei A and C; otherwise: 
in nuclei A and F) from the compound level as shown here 
(see text for details) . In each nucleus, a number of low-lying 
states with given spin J and parity vr is explicitely specified. 
Above the last state, transitions can be computed by inte- 
grating over nuclear level densities (shaded areas). In stellar 
cross sections a* all transitions are additionally weighted by a 
Boltzmann distribution factor depending on the stellar tem- 
perature, spin, and the excitation energy of the involved state 
(see Eq.El). 



II. SUPPRESSION OF THE STELLAR 
ENHANCEMENT 

A. Stellar reaction rates 

The stellar enhancement factor (SEF) / is defined as 
the ratio of the stellar rate r* relative to the ground state 
rate r^ ''- H 
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The rate r^'^^ derived from cross sections a^^^ measured 
in the laboratory is the same as r^'^' because so far all 
experiments use target nuclei in their ground states. The 
SEF is a measure of the influence of the thermally excited 
target states in the hot plasma. 

Astrophysical reaction rates are usually defined as giv- 
ing the number of a specific reaction occurring per time. 
Here, we constrain ourselves to two-body reactions of 
nuclei and nucleons. The concept of the stellar rates 
suppression introduced below is easily extended to other 
reaction types. Reaction cross sections are folded with 
the energy distribution of the interacting nuclei to ob- 
tain the reaction rate. The energy distributions of nuclei 



and nucleons in an astrophysical plasma follow Maxwell- 
Boltzmann distributions in most applications, thus yield- 
ing i 
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for reactions proceeding from target state i with reaction 
cross section ai, where iii, 712 are the number densities 
of the interacting nuclei, T is the plasma temperature, k 
denotes the Boltzma nn cons tant, and F is a renormal- 
ization factor F — y/8/{TTfi) with fi = A1A2/ {Ai + A2) 
being the reduced mass number A. 

When nuclei are in thermal equilibrium with their en- 
vironment, their excited states are populated according 
to a Boltzmann factor M 



{2Ji + l)e-T^ 
E„ (2J„ + l)e-i 



(3) 



with Pi, Ji, Ei denoting the relative population, spin, 
and excitation energy of state i, respectively. Each of 
the states is bombarded with Maxwell-Boltzmann dis- 
tributed projectiles which would require to have a sepa- 
rate rate for each target state weighted by the population 
factor of the state i from which the reaction proceeds. 
It was shown in [9] (see also [10]) that by making use of 
the reciprocity theorem for nuclear reactions and detailed 
balance (assuming thermalization of both initial and fi- 
nal states of a reaction), for compound reactions the rate 
equation can be simplified to 
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In order to avoid additional computations caused by the 
population coefficients and also to avoid having a temper- 
ature dependent stellar cross section, the effective cross 
section a'^^ was introduced above, which sums over all 
bound states in the initial and final system (denoted by 
1 and i , respectively; the energetics of the transitions is 
shown in Fig. [T]) 
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This is a theoretical construct (as any measurement 
would always proceed on a certain initial state and thus 
neglect the sum over target states) but it is useful in two 
respects. Firstly, it simplifies the computation of the rate 
and therefore is utilized in all astrophysical compound re- 
action codes. Secondly, it allows us to easily find a reci- 
procity relation between forward and inverse rate by re- 
membering that Ecr'^^ obeys reciprocity between forward 
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and inverse reaction due to detailed balance. It should be 
noted that only stellar reactivities TZ* (and thus stellar 
rates r*) obey reciprocity (as long as detailed balance is 
applicable) whereas rates derived from ground state cross 
sections a^^^ — J^j '^oj do not, unless the SEF is equal 
to unity in the given direction. 

For reactions 1 + 2 3 + 4 with target nucleus 1, 
projectile 2, final nucleus 3, and ejectile 4, the relation 
between backward and forward stellar reactivity is given 
by UM 



T?* — 



(2J0 + 1) (A1A2 



1) 



3/2 



— e -7^I 



(7) 



where Qn is the reaction Q value, J" denote ground state 
spins, and G are nuclear partition functions summing 
over states i and integrating over a level density p above 
the last discrete state m included [lOl. Ifll: 



G(r)= ^(2j, + i; 
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This partition function also appears in Eq. 2] where it is 
sufficient to compute it once and separately from the rate 
integration. 

Stellar capture reactions 1 + 2 — > 3 + 7 are related to 
stellar photodisintegration by [13, [HI 



7^;^ = (2j2"+i) 
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Reaction Q value and stellar enhancement 
factor 



Figure [T] shows a sketch of the energetically allowed 
transitions included in the effective cross section defined 
by Eq. ©. It is obvious that there are more transitions 
possible to and from states of the nucleus being the final 
nucleus in a reaction with positive Q value. Therefore, 
assuming a similar level structure in all involved nuclei, it 
is expected that the SEF (see Eq. [T]) of a given reaction 
will be smaller for the exothermic direction /forw than 
for the endothermic one /rev (here we define the forward 
reaction to be the one with positive Q value and the 
reverse reaction having negative Q value): 



/rev > /i 



for\ 
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This is especially pronounced in photodisintegration re- 
actions due to the many possible 7 transitions [l^, [ll] ■ 
In consequence, aiming at performing a measurement as 
close as possible to the stellar value, an exothermic reac- 
tion should be chosen. 



Another impact of the Q value is found by inspection 
of Eqs. ([7]) and ^ where the Q value appears in an expo- 
nential. For numerical consistency and to obtain proper 
equilibrium abundances when forward and reverse reac- 
tion are both fast and in equilibrium, reaction network 
codes avoid employing separate rates for the two direc- 
tions but rather make use of these equations. Taking an 
endothermic reaction as starting point for application of 
the equations would lead to a large value of the exponen- 
tial term, amplifying any numerical errors inherent in the 
original rate and in the Q value. This is mainly impor- 
tant when dealing with rate fits. In many astrophysical 
reaction network codes, the rates are implemented not 
as large tables but as fits with a smaller number of pa- 
rameters per reaction. Any deficiency in the fit would be 
amplified when computing an exothermic rate from an 
endothermic one. 

For the above reasons, it was commonly assumed that 
it is always preferable to determine the cross section and 
rate of an exothermic reaction and not those of an en- 
dothermic one. Here, we want to correct that notion by 
showing that there are cases for which 



J'rcv ^ /forw 
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The basic idea is to realize that although there are more 
transitions energetically possible to the final states of 
exothermic reactions, some of them may be suppressed 
and thus not contributing. Of course, it is obvious that 
not all transitions are of equal strength. Quantum me- 
chanical spin and parity selection rules and centrifugal 
barriers (or lack thereof) may prefer certain transitions 
over others. This will be important in reactions with 
small IQI and in nuclei with large level spacings. In both 
cases, only a small number of transitions will be possible 
and the spins can give larger weight to an even smaller 
subset. However, for reactions with sizeable Q values or 
involving nuclei with high level densities this suppression 
due to spins will not be sufficient because there will al- 
ways be a number of states with matching spins. 

Transitions to higher lying excited states of a nucleus 
proceed at lower relative energy. Except for s-wave neu- 
trons, also transitions at lower relative energy will be 
weaker than those at larger relative energy. If the sup- 
pression of transitions with smaller relative energy is dif- 
ferent in the entrance and exit channel of the reaction, 
this may also result in /rev < /forw The strongest sup- 
pression for charged particles is due to the Coulomb bar- 
rier. Having different Coulomb barriers in the entrance 
and exit channel, e.g. in (n,p) or (p,a) reactions, can 
more strongly suppress the transitions to the nucleus 
with higher Coulomb barrier than to the one with lower 
Coulomb barrier. With respect to Fig. [1] and assuming, 
e.g., a reaction A(n,Q;)F this means that most transitions 
to states in nucleus F are suppressed and the contribut- 
ing transitions may be fewer than those accessing states 
in nucleus A. 

This Coulomb suppression of the SEF is a general prin- 
ciple almost independent of nuclear structure and will act 
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FIG. 2: Reaction Q values for (p,n) and (a,n) reactions with 

frev ^ /forw ■ 



for a large range of nuclei. Whether the suppression is 
strong enough to yield /rev < /forw depends on the size 
of the Q value relative to the Coulomb barrier, i.e. the 
effect can occur in a reaction provided that there are dif- 
ferent Coulomb barriers in the entrance and exit channel 
and \Q\ is low compared to the Coulomb barrier. The 
strongest impact is to be expected when the forward reac- 
tion involves neutrons in the entrance channel which form 
a compound state by s- waves on excited target states and 
charged particles experiencing a high Coulomb barrier in 
the exit channel. As discussed in Sec. IIII Gl the reaction 
^^Rb(p,n)*^Sr is an excellent example for such a case. A 
quantitative exploration of the suppression effect across 
the nuclear chart is given in the following section. 

Not only theoretically interesting, this Coulomb sup- 
pression effect is also important for experiments because 
it allows to directly determine an astrophysically relevant 
rate by measuring in the direction of suppressed SEF. 
The above mentioned complication of fitting rates with 
negative Q values can be circumvented by directly apply- 
ing detailed balance and numerically computing the rate 
for the forward reaction before performing a fit. This is 
possible when /rev ~ 1- Subsequently, fits for both rates 
can be obtained in the standard way. As an example, an 
application of this procedure to the rate of ^^Rb(p,n)^^Sr 
is shown in Sec. IIII Gl 



C. Exploration of the SEF suppression across the 
nuclear chart 

In this section, we quantitatively study the SEF sup- 
pression introduced and discussed in the previous section. 
Using NON-SMOKER results we compared /forw 

and /rev for reactions involving light projectiles (nucle- 
ons, a) and targets from Ne to Bi between the proton 
and neutron driplines. To avoid trivial cases, only reac- 
tions with /forw//rev > 1-1 are considered. Furthermore, 
the SEF were computed for T < 4.5 GK to find cases 
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FIG. 3: (Color online) Targets for endothermic reactions with 
/rov < /forw in the nuclear chart, where charge is denoted by 
Z and neutron number by A'^. The reaction type is given by 
the label. Only capture or photodisintegration reactions are 
shown. Also printed for orientation are stable and longlived 
nuclides. 



important in most nucleosynthesis environments and to 
eliminate cases only occurring at high temperature. Be- 
cause of our aim to provide guidance for experiments, we 
further only focus on examples with /rev < 1-5. Even 
with these restrictions we find 1200 reactions exhibiting 
such a strong suppression effect that /rev < /forw 

To check the dependence on the Coulomb barrier. Fig. 
[5] shows the obtained range of Q values still yielding 
/rev < /forw as a function of target charge Z for (p,n) and 
(ck,n) reactions with negative Q values. It can be clearly 
seen that larger maximal \Q\ is allowed with increasing 
charge Z. The different increase in permitted maximal 
Ql is different for the two reactions, reflecting the differ- 
ence in the height of the acting Coulomb barriers. Below 
each maximally allowed |Q| for each given charge, there 
is a range of other values. This scatter is mainly caused 
by the available Q values (as defined by the masses of the 
nuclei) and not by other effects such as spins and parities 
of the involved nuclei. Although the strengths of the in- 
volved transitions also depend on spin and parity of the 
initial and final state, Coulomb repulsion dominates the 
suppression when the interaction energy is small, as it is 
the case for astrophysically relevant energies. 

Tables IB-fVl list the reactions found to have /rev < 
/forw according to the criteria discussed above. Stable or 
longlived target nuclei are marked specifically. Figures [3] 
and [5] locate the reactions in the nuclear chart. 

To current knowledge, not all shown reactions are of 
astrophysical importance. Others may be too far from 
stability to be accessible to experiments. Figure [3] dis- 
plays capture or photodisintegration reactions. Of par- 
ticular importance among them are the (a, 7) reactions, 
also given in Table HI They are located in a mass re- 
gion which is of interest in the p- or 7-process 0, Q. 
Although the 7 process synthesizes nuclei via photodis- 
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TABLE I: Targets for (a, 7) reactions 
but smaller SEP than their inverse 
longlived targets are in italics. 



with negative Q value 
reaction. Stable or 



3**Cd 


104 


130 Ce 


143 Eu 




153Lu 


l''3-y^ 


178p£ 




106 rpg 


131 Ce 


145 EU 


i5iEr 


155Lu 


1'4W 


179 p£ 


99ln 


109 rpg 


"2Ce 


"■^Eu 


i52Er 


156Lu 




189pt 


i°iln 


llOrpg 


133 Ce 


149 Eu 


153 Er 


162Lu 


i"W 


181pt 




"^Te 


134 Ce 


142 Gd 


i54Er 


163Lu 


161 Re 


182pt 


103 


112rpg 


135 Ce 


143 Gd 


i55Er 


166Lu 


i6''Re 


183 Pt 


105 


114rpg 




144 Gd 


156 Er 


i^'^Lu 


I'^iRe 


184 Pt 


107jj^ 


llSrpg 


i29pr 


145 Gd 




154Hf 


i'2Re 


185 Pt 


^^Sn 


iierpg 


i3ipr 


146 Gd 


158 Er 


156Hf 


i''5Re 


186pt 




"''Te 


i33pr 


"''Gd 


159 Er 


157jj£ 


i''9Re 


187p£ 


^^Sn 


113 J 


i35pr 


148 Gd 


i60Er 


158Hf 


164 Os 


"3Au 


^'"'Sn 


117j 


141 p J. 


149 Gd 


161 Er 


159jj£ 


166 Os 


"6Au 


lo^Sn 


119 J 


i42pr 






i6iHf 


I'^OOs 


182 Au 


i^^Sn 


lo^Xe 


132Nd 


'^^Gd 




162Hf 


1^2 Os 


184 Au 


i^^Sn 


"*Xe 


133 Nd 


i4^Tb 


153 Tm 


163 Hf 


l''30s 


185 Au 




"^Xe 


i34Nd 


148 Tb 


i55Tm 


164 jj£ 


1''40S 


189 Au 


lo^Sn 


i^'^Xe 


i35Nd 


149 Tb 


158 Tm 


165Hf 


l''5 0s 


ISOjjg 


I'^^Sn 


i2iXe 


i36Nd 


i42Dy 


i59Tm 


166 jj£ 


l''60s 


182 Hg 


lO'^Sn 


i^^Xe 


i37Nd 


144 Dy 


i6iTm 


167 jj£ 


I'^'Os 


184Hg 


I'^^Sn 


i23Xe 


13«Nd 


146Dy 


163 Tm 


169 jj£ 


l''80s 


185Hg 


i^^Sn 


i^iCs 


i39Nd 


"'Dy 


165 Tm 


i™Hf 


l''90s 


186 jjg 


"'^Sn 


Cs 




148Dy 


i52Yb 


171 JJ£ 


IBOQg 


187Hg 


"^Sn 


i25Cs 


i3ipm 


149Dy 


153 Yb 


172 jj£ 


18l0s 


188 Hg 




i27Cs 


133 Pm 


i50Dy 


154 Yb 


173 jj£ 


1820s 


189 jjg 




i23Ba 


i37pm 


151 Dy 


156 Yb 


155 Ta 


1830s 


196 Hg 


"^Sn 


i^^Ba 


149 


i52Dy 


157Yb 


165 Ta 




i9iHg 




i25Ba 


i34Sm 


154 Dy 


158 Yb 


i6''Ta 


172 


192Hg 


Sb 


i26Ba 


i36Sm 




159Yb 


169 Ta 


175 jj. 


193 Jjg 


lO'^Sb 


i^'^Ba 


"^Sm 




160Yb 


158-yy 


177 


194 Jjg 




i^^Ba 


i3*Sm 




161 Yb 


160W 


179 jj. 


195 jj 


io»Sb 


i29Ba 


i39Sm 




162Yb 


164W 


i8ilr 




lo^Sb 




"OSm 


149Ho 


163 Yb 


166W 


183 Ir 




"iSb 




"iSm 


150Ho 


164 Yb 


i6rw 


i8'lr 


•'''"Hg 


"^Sb 


i29La 


i42Sm 


i5iHo 


165 Yb 


168W 


168pt 


''"''Hg 


"3sb 


i3iLa 


143 Sm 


154Ho 


166 Yb 


169-VV 


170p^ 


187rj,j 


"5sb 


Ce 




155HO 


167Yb 




l'4pt 


189rjiJ 


"3Sb 


128 Ce 




i5'Ho 




I'l-W 


l'6pt 


192rj,j 


102 Te 


129 Ce 


i3^Eu 


159Ho 


178 Yb 


1'2W 


l^'Pt 





TABLE 11: Targets for (p,7) 
but smaller SEP than their 
longlived targets were found. 



reactions with negative Q value 
inverse reaction. No stable or 



18 Ne 


39y 


53Ni 


^6Sr 


i03Sn 


ii4Xe 


152Yb 


180Hg 


24Si 


44Cr 


54Ni 


86Ru 


i04Sn 


i26Nd 


153 Yb 


^^^Hg 


29s 


43 Mn 


51 CU 


88Ru 


i04Te 


i30Sm 


154 jj£ 


186pb 


33 Ar 


46Fe 


56 Zn 


90pd 


106 Te 


136Gd 


157jj£ 


188pb 


31 K 


4^Fe 


5'Zn 


92pd 


108 Te 


138Dy 


158w 




36Ca 


48Fe 


58Zn 


95 Cd 


109 Te 


144Dy 


160w 




3'Ca 


49Fe 


69Ge 


96 Gd 


iiOTe 


i48Er 


1640S 




38Ca 


4'Co 


61 Ge 


10° Sn 


ii2Xe 


isoYb 


i™Os 




46Ti 


52Ni 


62Ge 


102 Sn 


ii3Xe 


isiYb 


"6pt 
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PIG. 4: (Color online) Targets for endothermic reactions with 
/rev < /forw in the nuclear chart, where charge is denoted by 
Z and neutron number by A''. The reaction type is given by 
the label. Only reactions without 7 channels are shown. Also 
printed for orientation are stable and longlived nuclides. 



TABLE HI: Targets for (7,n) reactions with negative Q value 
but smaller SEP than their inverse reaction. No stable or 
longlived targets were found. 



36Ne 


51p 


™Ca 


102 Zn 


125y 


i3^Nb 


186 Ce 


188 Nd 


34Na 


53 p 


82 Cr 


ii6Se 


131y 


130Mo 


i88Ce 


194Nd 


42Mg 


55 p 


83 Mn 


i2iRb 


i26Zr 


140Mo 


i8'Pr 


236 jjf 


40 AI 


58 g 


90pe 


i22Sr 


1302, 


i5iRh 


i89pr 


2460s 


50Si 


64 Ar 


94Ni 


i28Sr 


129Nb 


152pd 


I93pr 


1^6 Au 



integrations, it becomes obvious that even for a captures 
with negative Q value the SEF of the capture reaction is 
smaller than the one of the photodisintegration. Thus, 
a measurement of the capture includes more astrophysi- 
cally relevant transitions and is closer to the stellar value 
than a measurement of the photodisintegration in the 
laboratory. Many of the cases with negative Q value are 
found at stability, providing interesting examples for ex- 
perimental study. A second class of interesting reactions 
are the (p,7) and (7,p) reactions in Fig. [3] and Tables [TTl 
IIVI They involve unstable targets and are of interest in 
the rp process jll] and the vp process [l6| . 

Figure |4] displays reactions without a photon channel. 



TABLE IV: Targets for (7,p), (n,a), and {n,7) reactions with 
negative Q value but smaller SEF than their inverse reaction. 
No stable or longlived targets were found. 



(7,P): 
(n,a): 

(11,7): 



28 Ar 
45 AI 



4''Mn 
62 Ar 



83 



Cr 



132 



Pd 



=Au 



"Cd 
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TABLE V; Targets for (p,ii) reactions with negative Q value 
but smaller SEP than their inverse reaction. Stable or 
longlived targets are in italics. 



^^Si 


''As 


"'iRh 


""Sb 


144 


""Tb 


i^^Lu 


'■"Pt 


«Ar 




'"■'Rh 


125Sb 


'^'^Nd 


161 




198p^ 


41 K 


''Br 


Rh 


128 j,^ 


148^^ 






200pt 








130 j,^ 


""Nd 


166Dy 


182jjf 


195 






""Pd 


132 


'^'Pm 


163Ho 


1^8 Ta 


'^'Au 


4^Sc 




""Pd 


127-^ 


'^'Pm 






«8Au 








129-|- 


"^Sm 








Sly 


»^Rb 






"*Sm 




185W 


''"'Hg 




*Sr 


">"Ag 




156 Sm 




186 W 


202jjg 


eope 




"^Cd 




"^Eu 


le^Tm 


188 W 


''"^Hg 


/l/ TIT ■ 


Q/ rr 


'"'Cd 


"^Cs 


1 '^1 r-i 

'•"Eu 


'"■^Tm 


"'•'Re 


900 ml 


66 Ni 


""Zr 


"«Cd 






1^1 Tm 


'^'Re 






■''Nb 


'"In 


CS 


155EU 


"'Yb 


'""Os 




e^Cu 


'"Mo 


"'In 






112 




205 


™Zn 


'""Mo 


'''"Sn 


"^La 




myb 


"*Os 




^^Zn 


99 j,^ 


''■'■'Sn 




'""Gd 


176 Yb 


189jj. 




^'Ga 




'''''' Sn 




153Tb 


178Yb 


191 




^'Ge 


'"^Ru 


i26Sn 


"^Ce 


155Tb 


"3Lu 


193 jj. 




■''As 




''''Sb 


141 


157Tb 




196p^ 





TABLE VI: Targets for (n,p) reactions with negative Q value 
but smaller SEP than their inverse reaction. No stable or 
longlived targets were found. 





«3Ar 


102 Zn 


""Zr 


153pd 


216Dy 


243^ 


250Os 


46gi 


««Ca 


121 Kr 


"OMo 


i»5Ba 


224Er 


2«Re 


269pi3 


50Si 


™Ca 


i^^Sr 


i52pd 


i9«Nd 


230Yb 


2480s 





Mostly interesting is the large number of (p,n) reactions 
located along stability (see Table W\ . A recent investi- 
gation of the 7 process has pointed out the importance 
of (n,p) reactions close to stability According to our 
findings, it is best to study the endothermic (p,n) reac- 
tion when trying to experimentally constrain the astro- 
physical rate. Our measurement of ^^Rb(p,n)®^Sr, de- 
scribed in the following section, is an example for such 
an experiment. 

Of minor or no astrophysical relevance are the (Q;,n) 
and (n,p) reactions in Fig. [Hand Tables IVlIl IVII because 
they involve neutron-rich targets and these reactions are 
much slower than other possible reactions on the same 
nuclei. 

Figure [3] and Tables EH |TV] list (n,7) and (7,n) reac- 
tions on very neutron-rich nuclei far off stability. They 
are relevant in the r process [13, [3- Obviously, the 
suppression of the SEF is not caused by the Coulomb 
barrier. The nuclear structure (spins and parities of ex- 
cited states and the nuclear level density) is important 
in those cases, as the level density is low and transitions 
are favored or suppressed by selection rules and centrifu- 
gal barriers. Contrary to the Coulomb suppression of 
the SEF, this type of suppression is strongly dependent 
on the nuclear spectroscopy assumed in the calculation. 



TABLE VII: Targets for (Q,n) reactions with negative Q 
value but smaller SEP than their inverse reaction. Stable 



or longlived targets are in italics. 







IOSRu 


133 Te 


i™Pm 


"5Tm 


190W 


198p^ 


32Si 


>"'Kr 


108 Ru 


i3*Te 


1153?^ 


Tm 


192-vv 


199p^ 


35 p 




hORu 


135 Te 


ii^-ipm 


"8Tm 


197w 


2"Pt 


36 g 


»^Kr 


ii^Ru 


136 Te 


ii^^Pm 


"^Tm 


200^ 


2°lpt 




»»Kr 


ii^Rii 


i3^To 


'■'''Sm 


i8»Tm 


201-w 


202pt 


''Ar 


««Kr 


""Rh 


138 Te 


i56Sm 


i8iTm 


202^ 


201pt 


'">Ar 




ii-'Rh 


129j 


i58Sm 


i82Tm 


203^ 


205p^ 


^"Ar 




i"=Rii 


131 J 


155Eu 


183 Tm 


206^ 


206pt 


40k 


^■iRb 


lOBRh 


132j 


i5^Eu 


i8iTm 


207^ 


210pt 


41k 


x''m 


lO^Rh 


133j 


158Eu 


185 Tm 


210w 


212pt 


43 k 


»«Rb 


lOSRh 


135j 


i5fEu 


i86Tm 


212w 


213pt 


^'■Ca 


»^Rb 


iiiRh 


136j 


i«°Eu 


i88Tm 


211W 


218pt 




****Rb 


i"Rh 


137j 


i«2eu 


isiTm 


l™Ro 


220pt 




^ORb 






163Eu 


135 Tm 


i3iRo 


221pt 


^'Sc 


^iRb 


'""Pd 




I'^iEu 


i"Yb 


132Ro 


222pt 




^^Sr 


""Pd 




1'35Eu 


178Yb 


i9iRc 


223pt 


«Sc 


«=Sr 


iiipd 


i38Xc 


166Eu 


i86Yb 


200RC 


221pt 




""Sr 


ii^Pd 


i39Xc 


le^Eu 


188 Yb 


201 Rc 


131 Au 






ii-ipd 


"OXc 


i«»Eu 


190Yb 


209 Rc 


188 Au 


62 


™Sr 


ii«Pd 


"iXo 


i™Eu 


19lYb 


2iORe 


is'^'Au 


51y 


31 Sr 


'"•'Ag 


"=Xo 


1^3Eu 


194Yb 


2iiRe 


201 Au 


53y 


3=Sr 


iiiAg 


i«Xo 


"'Gd 


196Yb 


212Ro 


205 Au 


=2Cr 




ii^Ag 




""Gd 


198Yb 


2i3Re 


213Au 






ii=Ag 


"^Cs 


i«iGd 


200Yb 


215Ro 


2iiAu 


5«Cr 


"BSr 


""Ag 


i38Cs 


le^Gd 


"8Lu 


2"Ro 


215Au 


•'^ Mn 


86y 


ii^Ag 


13fCs 


i^iQd 


"9Lu 


18l0s 


216Au 




OOy 


'"Cd 


"°Cs 


i^^Tb 


183Lu 


""Os 


2"Au 


^^Mn 


Sly 


"^Cd 


"iCs 


163:pjj 


i8iLu 


l»30s 


218Au 




92y 


""Cd 


i«Cs 


if^iTb 


185Lu 


i^iQs 


218Au 




93y 


ii»Cd 


'"'Ba 


166:pjj 


186Lu 


l«50s 


220 Au 


«iCo 


94y 


i20Cd 


"'Ba 


les^b 


i8''Lu 


l«60s 


221 Au 




95y 


i"ln 


i-'OBa 


if^^Tb 


188Lu 


l«80s 


222 Au 




97y 


ii^In 


i-iiBa 


170 Tb 


189Lu 


i««Os 


223 Au 


66 Ni 


""Zr 


i^iln 


1*2 Ba 


i^iTb 


i9iLu 


200Os 


221 Au 


68 Ni 


'"'Zr 


i^^In 


"3Ba 


173Tb 


135Lu 


2«20s 


225 Au 


«^Cu 


'"'Zr 


i==In 


"''Ba 


175Tb 


i3^Lu 


2«30s 


226 Au 


6«Cu 


""Zr 


i^'ln 


"=Ba 


176Tb 


138Lu 


2«l0s 


22^Au 




""Zr 


""Sn 


i«Ba 


179Tb 


2»iLu 


2«80s 


228 Au 


™Zn 






i"Ba 




2»3Lu 


209Os 


220 Au 


^^Zn 






"8Ba 


166Dy 


183Hf 


210Os 


23lAu 




1002, 


i2=Sn 


i™Ba 


i"Dy 


i8iHf 


21l0s 


2iiAu 


■"Ga 


1022, 


i^i^Sn 


1-13 La 


i^^Dy 


i™Hf 


2130s 


""Hg 


"Ga 


"■iNb 


i^'-Sn 


i-i-iLa 


i8iDy 


192gf 


21l0s 


'■"'■'Hg 


^^Ga 


"^Nb 


i^^Sn 


1-18 La 




i35Hf 


2150s 


'■"'*Hg 


''^Ge 


Sf^Nb 


i^^Sn 


i™La 


I'^^Ho 


196 


21«Os 


205 Hg 


™Ge 


s^Nb 


i3"Sn 


i^^Ce 


i™Ho 


197jjf 




206Hg 


™Ge 


38 Nb 


i^iSn 


118 Ce 


"2 Ho 


198jjf 


""Ir 


207Hg 


^»Go 


™Nb 


i^^Sn 


i™Ce 


"3 Ho 


199jjf 


194 J, 


208Hg 


"As 


lo^Nb 


i^^Sn 


152 Co 


"iHo 


200jjf 


195j, 


200Hg 




lo^Nb 


i^iSn 


151 Ce 


"5 Ho 


201jjf 


196 


210Hg 


™As 




i^f'Sn 


ii^Pr 


"''Ho 


204 jjf 


197j^ 


212Hg 


**iAs 


""Uo 


i^^Sn 


118 Pr 


i"Ho 


i85Ta 


199 J, 


""Hg 


™Se 


"'"Mo 




I58pr 


"8ho 


i88Ta 


203j, 


215Hg 


™Se 


i«iMo 


i2=Sb 


"'"Nd 


i™Ho 


i89Ta 


2iilr 


216Hg 




i°2Mo 


i^'Sb 


i5iNd 


i8iHo 


i9«Ta 


2i3lr 


222Hg 


Se 


Mo 


Sb 


Nd 


1S2tt 

Ho 


Ta 


214t 

Ir 


224 TT 

Hg 




i™Mo 


i^^Sb 


i5iNd 


183Ho 


i^^Ta 


2i5lr 


226Hg 




i™Mo 


i^«Sb 


""Prn 


i8iHo 


i°8Ta 


216j^ 


229 Hg 


"'Br 


'""Tc 


i^iSb 


152 Pm 


185Ho 


i^^Ta 


2"lr 


23lHg 


''Br 


101 j,^ 


i^^Sb 


153 Pm 


™Er 


™iTa 


2i8lr 


233 jj^ 




lo^Tc 


i^^Sb 


I5ipm 


"^Er 


2"iTa 


219j, 




83b^ 


103 


"'■<'' Te 


I55pm 


i8iEr 


205 Ta 


220j, 






105-pc 




I56pm 


i85Er 


20^Ta 


22ilr 








130 j,^ 


i5'Pm 


i86Er 


2°3Ta 


223j, 




''Kr 


">iRu 


i^iTe 


I58pm 


i^iEr 


188W 


225 J, 




'■'Kr 


i°^Ru 


i32to 


I59pm 


i™Er 


189W 


W6p^ 





TABLE VIII: Targets for (ci,p) reactions with negative Q 
value but smaller SEF than their inverse reaction. Stable 
or longlived targets are in italics. 



85Rb(p,n)85Sr 



39 Ca 


^''Ga 


^^Kr 


8«Zr 


97pd 


loiSb 


iiiCs 


134 Gd 




s^Ge 




87 Zr 


97 Ag 


i05sb 


120 Cs 


139Tb 




'^"Ge 


^«Kr 


84Mo 


98 Ag 


lo^Sb 


ii^Ce 


138Dy 




65 Qg 


™Rb 


86 Mo 


98 Cd 


108 Sb 


ii8Ce 


i45Tm 


58 CU 


««Ge 


^«Sr 


88 Mo 


199Cd 


i09sb 


i20Ce 


i46Tm 


«»Cu 


™As 


89Sr 


89 Mo 


103 Cd 


191 Te 


i25pr 


i47Tm 


58 Zn 


«8Se 


81 Sr 


9iMo 


99ln 


192 Te 


i27pr 




59 Zn 


69 Se 


83Sr 


92Ru 


i93ln 


194 Te 


i24Nd 




"OZn 


^OSe 


89 Zr 


93Ru 


i«In 


196 Te 


126Nd 




eiZn 


"Se 


82 Zr 


94Ru 


97Sn 


iiij 


i29pm 




«2Zn 


^^Se 


83 Zr 


95pd 


99Sn 


i96Xe 


i39Sm 




«Ga 


^^Br 


84 Zr 


96pd 


Sn 


ii^Xe 


135Eu 





Moreover, the Hauser-Feshbach model of compound re- 
actions may not be applicable anymore for the nuclei in- 
volved [l^ . Additionally, the individual rates are not im- 
portant in r process studies producing highly unstable, 
neutron-rich nuclei close to the dripline in an equilibrium 
between capture and photodisintegration [13, [2l| . 
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100 % 



85Srm 
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FIG. 5: Simplified decay scheme of the products of the 
85Rb(p,n)85Sr reaction. The half-lives of the reaction prod- 
ucts, the branching ratios, the spin and parity of the levels 
and the transitions used to determine the reaction cross sec- 
tion are indicated [231. 



III. EXPERIMENTAL STUDY OF 85Rb(p,n)85Sr 
A. General 

As an example of the suppression effect and for the 
derivation of the astrophysical rates for an endother- 
mic reaction, we experimentally studied the reaction 
85Rb(p,n)85Sr. 

Reactions of the (n,p) type have been shown to be im- 
portant in the 7 process [6| . This nucleosynthesis process 
creates proton-rich isotopes of elements beyond Fe which 
are not made in the s and r processes. It was shown 
to occur in hot 0/Ne layers of massive stars, either in 
a core collapse supernova explosion when the shockfront 
is passing these layers or already pre-explosively depend- 
ing on the initial mass of the star . At temperatures 
T > 2 GK photodisintegrations can act even within the 
short timescale of an explosion. The reaction sequences 
initially drive material from the bottom of the valley 
of stability to the proton-rich side by (7,n) reactions. 
Charged-particle emitting (7, a) and (7,p) reactions can 
deflect the reaction path to lower charge number. The- 
oretical investigations show that (7,n)/(7,p) branchings 
play a key role in the production of the lighter p nuclei 
whereas (7 ji)/ (7, a) branchings are important at higher 
masses [a, Q. Further reactions with the emitted neu- 
trons are mainly important in the freeze-out phase when 
photodisintegration ceases . The flow back to sta- 

bility is sped up by (n,p) reactions which are faster than 
f3 decays close to stability 0. Even at stability, (n,p) re- 
actions act to push material to lower proton numbers. In 
this context, *^Rb(p,n)*^Sr is directly important because 



it is the inverse reaction to ^^Sr(n,p) and we found that 
its SEF is smaller than the one of its inverse, despite of 
its negative Q value (see Sec. IIIIGp . 

The reaction ^^Rb(p,n)*^Sr is also important to test 
the predictions of astrophysical rates and their under- 
lying nuclear properties. Although many (n,p) and 
y) reactions important in the 7, rp, and lyp processes 
occur far from stability, the models and as- 
sumptions used in the prediction of the rates can be 
checked at stability. Especially suited for testing the re- 
liability of the optical potential used for the calculation 
of transitions involving protons are (n,p) and (p,n) reac- 
tions because the proton width is smaller than the neu- 
tron width at practically all energies (except very close 
to the neutron threshold) and thus determines the cross 
section. 

We measured *^Rb(p,n)^'^Sr using the activation 
method. Thin RbCl targets were bombarded by proton 
beam provided by the Van de Graaff and cyclotron accel- 
erators of ATOMKI 0. The (p,n) reaction on ^^Rb can 
populate both the ground and metastable states of ^^Sr 
[23 |. To determine the cross section of the ®^Rb(p,n)®^Sr^ 
reaction the 514.01 keV 7 line was used, in the case of 
the ^^Rb(p,n)*^Sr™ reaction cross section the yield of the 
231.84 keV transition was measured. In the following 
Sees, nil B1 nil Dl a detailed description of the experimen- 
tal technique is given, while the experimental results are 
given in Sec. IIIIEI A comparison to theory and the final 
astrophysical reaction rates are provided in Sees. IIIIFI 
and lmUl 
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FIG. 6: Typical 7-spectra taken after the irradiation of RbCl targets with 2.4 (left panel) and 3.8 MeV (right panel) proton 
beams. The 514.01 keV peak from the *^Rb(p,n)*^Sr® reaction can be well separated from the annihilation peak as can be 
seen in the insets. The length of the waiting time (t^,) between the end of the irradiation and the start of the 7-countings 
were 540 (Ep = 2.4 MeV) and 30 min (Ep = 3.8 MeV). The lower panels show typical spectra taken in the repeated activity 
measurement approximately one month after the irradiations (for details see text). 



B. Target properties and the determination of the 
number of target atoms 

The targets were made by evaporating chemically pure 
(99.99%) RbCl onto two different kinds of Aluminum 
foils: the thicker one had a chemical purity of 99.999% 
and thickness of 50 fxin, while the purity and the thick- 
ness of the thinner one was 99% and 2.4 /zm, respectively. 
The distance between the evaporation boat and the tar- 
get backing was 10 cm, therefore it was possible to as- 
sume that the evaporated layer is homogeneous. This 
assumption was proved using Rutherford Backscattering 
Spectroscopy (RBS, see later). Targets with different 
thicknesses were used, thicker ones (on thicker backings) 
were employed for irradiations at lower and thinner ones 
(on thinner backing) at higher energy. Owning to this 
treatment, the yield of the investigated 514.01 keV peak 
was always higher than, or comparable to that of the 511 
keV annihilation peak - and this way the separation of 
the peaks was achieved - as it is demonstrated in the 
upper part of Fig. [S] 

The number of the target atoms was determined with 
Rutherford Backscattering Spectrometry (RBS) at the 
Nuclear Microprobe facility of ATOMKI (^,111, 113 • As 
a consistency check, in the case of the targets evaporated 
onto the thinner backing, weighing was also used to de- 
termine the number of target atoms. The weight of the 
Al foil used as backing was measured before and after the 



evaporation and from the difference — assuming that our 
target is uniform — the number of target atoms was cal- 
culated. The results of the two different methods used to 
determine the number of target atoms are in very good 
agreement (< 3% difference). 



C. Irradiation 

The RbCl targets were bombarded with a proton beam 
provided by the Van de Graaff and cyclotron accelerators 
of ATOMKL The energy of the proton beam was between 
2 and 4 MeV, this energy range was covered in 200 keV 
steps. The beam current was typically 600 nA. Each 
irradiation lasted approximately 7 — 8 hours. The low 
energy irradiations (2.2 MeV, 2.4 MeV and 2.6 MeV) 
have been carried out using the Van de Graaff. At and 
above 2.6 MeV the cyclotron accelerator was used. The 
cross section at Ep — 2.6 MeV was measured with both 
accelerators and no difference was found. 

An ion implanted Si detector was built into the irradi- 
ation chamber at 6* = 150° relative to the beam direction 
to measure the yield of the backscattered protons. The 
backscattering spectra were taken continuously and were 
used to monitor the target stability. Having a beam re- 
stricted to 600 nA, no target deterioration was found. 

For calculating the reaction cross section the proper 
knowledge of the incident particle flux is necessary. To 
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obtain this, the collected charge was measured in a cham- 
ber similar to the one in After the beam defining 
aperture, the whole chamber served as Faraday cup to 
collect the accumulated charge. A secondary electron 
suppression voltage of —300 V was applied at the en- 
trance of the chamber. The beam current was kept as 
stable as possible but to follow the changes the current 
integrator counts were recorded in multichannel scal- 
ing mode, stepping the channel in every minute. This 
recorded current integrator spectrum was used for the 
analysis solving the differential equation of the popula- 
tion and decay of the reaction products numerically. 



D. Activity determination 

Figure [5] shows the simplified decay scheme of the 
85gj-g,m jgQtopes. To determine the cross section of 
the ®5Rb(p,n)®5Srf reaction the 514.01 keV, for the 
^^Rb(p,n)^^Sr™ reaction the 231.84 keV gamma hue was 
used. 

For measuring the induced 7-activity a lead shielded 
HPGe detector was used as in our previous (p,n)-study 
[ij. After each irradiation, a cooling time of one hour was 
inserted in order to let the disturbing shortlived activities 
decay. The 7 spectra were taken for 12 hours and stored 
regularly in order to follow the decay of the short-lived 
reaction product. 

Figure [6] shows typical spectra collected after irradi- 
ating RbCl targets with the 2.4 MeV (left panel) and 
the 3.8 MeV (right panel) proton beams. The yield of 
the 511 keV peak was always less than or comparable to 
the investigated 514.01 keV transition, as shown in the 
insets. The ^^^Si has a relatively long halflife (T1/2 = 
64.84 d). Due to this, the activity measurement could be 
repeated for each target after approximately one month, 
when the intensity of the 511 keV radiation was substan- 
tially reduced. The spectra taken in the repeated ac- 
tivity measurement for the 2.4 and 3.8 MeV irradiations 
are shown in the lower panels of Fig. [Sj The two mea- 
surements yielded consistent cross sections proving the 
proper separation of the 511 keV and 514.01 keV peaks. 



E. Experimental results and comparison with 
literature data 

In the case of the *''^Rb(p,n)^^Sr^ reaction, two sepa- 
rated analysis were done. The agreement between the 
cross sections derived in the 7-counting after the irradi- 
ation and the ones from the repeated activity measure- 
ment was always better than 4%. The final results were 
calculated from the average weighted by the statistical 
uncertainty of the two 7 countings. The halflife of the 
^^Sr"* is shorter, therefore the yield of the 231.64 keV 7 
radiaton was measured only after the irradiation. The fi- 
nal experimental result can be found in Table HXl Partial 
cross sections leading to the ground and isomeric state of 
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FIG. 7: Experimental (full triangles) and theoretical (lines) 
astrophysical S factors of *^Rb(p,n)*^Sr. The solid line is the 
S factor calculated with the modified proton optical potential 
introduced in pj and the dashed line shows the result using 
the standard proton optical potential from [3l|] with low en- 
ergy modifications by [3^ (see text). 

^^Sr can be found in The error of the cross section 
values is the quadratic sum of the following partial errors: 
efficiency of the HPGE detector system (6%), number of 
target atoms (< 3.3%), current measurement (3%), un- 
certainty of the level parameters found in the literature 
(< 4%) and counting statistics (< 4%). The quoted er- 
rors of the energies include the energy loss in the targets 
calculated with the SRIM code [13] , as well as the energy 
stability of the cyclotron and Van de Graaff accelerators. 

The measured total cross sections cover 3 orders of 
magnitude, varying from 0.06 to 20 mb. Table HXl lists 
the measured cross sections a and the S factors, the latter 
being defined as [3] 

S{E)^^e-^-^ , (12) 

with the Sommerfeld parameter 77 accounting for the 
Coulomb barrier penetration. 

The cross section of the *^Rb(p,n)*^Sr reaction was al- 
ready investigated by [l^ between Ec.,„. = 3.1 and 70.6 
MeV. However, their accuracy is not sufficient for astro- 
physical applications, mainly because of the large uncer- 
tainty of the cm. energies. Moreover, there is only one 
data point in the relevant energy region for the 7 process 
and it bears an uncertainty of ± 0.5 MeV in the cm. 
energy. 



F. Comparison to theory and implications for the 
proton optical potential 

The measured S factors are compared to theoreti- 
cal predictions obtained with the code NON-SMOKER 
[ll].ll^ in Fig.[71 The standard calculation applied a pro- 
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TABLE IX: Details of the irradiations and the resulted cross sections (astrophysical 5* factors). 





Ec.m. 


Accelerator 


Collected 


Total <T 


S factor 


[MeV] 


[MeV] 




charge [mC] 


[mharn] 


[iQ^MeVharn] 


2.20 


2.16 ± 0.008 


Van de Graaff 


48.27 


0.058 ± 0.006 


7.13 ± 0.67 


2.40 


2.34 ± 0.016 


Van de Graaff 


37.23 


0.224 ± 0.019 


11.22 ± 0.96 


2.60 


2.57 ± 0.026 


Van de Graaff 


33.74 


0.582 ± 0.055 


11.01 ± 1.04 


2.58 


2.55 ± 0.027 


Cyclotron 


35.44 


0.569 ± 0.051 


11.35 ± 1.02 


2.79 


2.77 ± 0.028 


Cyclotron 


30.83 


1.20 ± 0.12 


10.65 ± 1.05 


2.98 


2.96 ± 0.030 


Cyclotron 


23.66 


2.12 ± 0.21 


9.59 ± 0.92 


3.18 


3.16 ± 0.032 


Cyclotron 


23.43 


3.77 ± 0.35 


9.40 ± 0.87 


3.37 


3.34 ± 0.035 


Cyclotron 


21.45 


5.66 ± 0.54 


8.37 ± 0.79 


3.57 


3.55 ± 0.036 


Cyclotron 


20.64 


9.60 ± 0.87 


8.32 ± 0.76 


3.76 


3.75 ± 0.037 


Cyclotron 


17.16 


14.31 ± 1.22 


7.73 ± 0.66 


3.96 


3.95 ± 0.040 


Cyclotron 


11.51 


19.65 ± 1.82 


6.93 ± 0.64 



ton optical potential widely used in astrophysical applica- 
tions, based on a microscopic approach utilizing a local 
density approximation [sH. Low-energy modifications, 
which are relevant in astrophysics, have been provided 
by 13^. As can be seen in Fig. [3 the theoretical energy 
dependence of the resulting S factor is slightly steeper 
than the data, although there is general agreement in 
magnitude. In the energy range covered by the measure- 
ment, the proton width is smaller than the neutron width 
(except close to the threshold) and thus uncertainties in 
the description of the proton width (and proton transmis- 
sion coefficient) will fully impact the resulting S factor. 
A recent investigation [ij] suggested that the strength of 
the imaginary part of the microscopic potential should be 
increased by 70%. We find that the energy-dependence 
of the theoretical S factor is changed in such a way as 
to show perfect agreement with the new data, as seen in 
Fig. [71 This independently confirms the conclusions of 
previous work pj. 



G. Astrophysical reaction rates 

Regarding the Coulomb suppression effect, a compar- 
ison of 1.03 < /pn < 1.08 and 2.6 < /np < 3.9 shows 
that the transitions to excited states of ^^Sr are more 
important than those to states in ^^Rb in the relevant 
plasma temperature range of 2 < T < 4 GK. The almost 
negligible stellar enhancement /p,i is due to the suppres- 
sion of the proton transmission coefficients to and from 
the excited states of ^'^Rb for small relative proton ener- 
gies because of the Coulomb barrier. There are only few 
transitions able to contribute due to the low Q value. As 
shown by the small /pn, the transition from the ground 
state of *^Rb dominates the proton channel. Obviously, a 
Coulomb suppression is not present in the neutron chan- 
nel. On the contrary, for this reaction /np is even more 
enhanced due to the spin structure of the available nu- 
clear levels and especially the large spin of *^Sr^. Be- 
cause of its large spin, it is connected to the (dominat- 
ing) low spin states in ®^Rb through higher partial waves 
than the excited states, such as the isomeric state, which 



have lower spins. Thus, the transitions from the ground 
state are suppressed by the centrifugal barrier relative to 
transitions from excited states and the latter will quickly 
become important, even at low temperature. As a conse- 
quence of the enhancement of /np and the suppression of 
/pn, it is more advantageous to measure the (p,n) direc- 
tion. Important transitions to states in ®^Sr are included 
in our data and the small impact of transitions from ex- 
cited states in *^Rb is within the experimental error. 

Applying Eq. [2] directly with the experimental cross 
sections already yields the stellar rate because the SEF 
is small in the (p,n) direction. The stellar rate of the 
exothermic (n,p) reaction can then be computed using 
Eq. ([7]). By computing the forward rates directly from 
the backward rates without using fits, the complication 
with the negative Q value in fitted data is also avoided. 

Table |X] gives the stellar reactivities (as defined by 
Eq. HI) for 85Rb(p,n)85Sr as well as for 85Sr(n,p)8^Rb. 
Our data covers an energy range sufficient to compute 
the rates between 2 and 4 GK. Because of the excellent 
agreement of theory with experiment (using the newly 
modified potential of [l|), we supplement the data with 
the theoretical values to compute the reactivities at lower 
and higher temperatures, applying the same errors as for 
the data. 

It is to be noted that fits of the rates should be ob- 
tained by first fitting the (n,p) rate and then deriving 
the (p,n) rate fit by modifying the fit coefficients accord- 
ing to detailed balance as given in Eq. (0) (see Til for 
details). For convenience, we provide the fit coefficients 
(including a 10% error) for the (n,p) reactivity in the 
widely used REACLIB format 



NpJZ* — exp oo + — 



Oi 



0-2 



1/3 



1/3 



04 Tg 



-a^T^^^ + aQ InTg 



(13) 



where JVa is Avogadro's number and the plasma tem- 
perature Tg — T/10^, with T in K. Using the usual di- 
mension of cm'^ s^^ mole^^ for Nj^TZ* the fitted coeffi- 



cients evaluate to oq 
as = 40.7296, 03 = 



33.2271+j"o:g, fli 
-67.9553, 04 = 



-0.886129, 
6.54471, as = 
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TABLE X: Astrophysical reactivities NaTZ* of the reactions 
**^Rb(p,n)**^Sr and **^Sr(n,p)**^Rb computed from experimen- 
tal data. The values in italics are at temperatures where the 
experimental data mostly contribute to the rate. The other 
values are computed by supplementing theoretical cross sec- 
tions using the modified optical potential. 



Temperature **^Rb(p,n)*^Sr **^Sr(n,p)**^Rb 
[10'' K] [cm^s'Vole'^] [cm'^s'^molc'^] 



0.10 


(1.72± 


0.17) xlO"'''' 


(1.19± 


0.2) X 10* 


0.15 


(2.21± 


0.22) X 10"^* 


(1.49± 


0.15) X 10* 


0.20 


(8.33± 


0.83) X 10"*^ 


(1.74± 


0.17)xl0* 


0.30 


(3.36± 


0.33) X 10"^^ 


(2.15± 


0.22) X 10* 


0.40 


(2.26± 


0.23) X 10'" 


(2.55± 


0.26) X 10* 


0.50 


(1.18± 


0.12)xl0~^'* 


(2.99± 


0.30) X 10* 


0.60 


(1.74± 


0.17)xl0~" 


(3.49± 


0.35) X 10* 


0.70 


(3.31± 


0.33) xlO"'' 


(4.09± 


0.41) X 10* 


0.80 


(1.77± 


0.18)xl0-^ 


(4.80± 


0.48) X 10* 


0.90 


(4.04± 


0.40) X 10"® 


(5.62± 


0.56) X 10* 


1.00 


(5.07± 


0.51)xl0~^ 


(6.57± 


0.66) X 10* 


1.50 


(1.28± 


0.13)xl0"^ 


(1.35± 


0.14) X 10^ 


2.00 


(8.30± 


0.83) 


(2.56± 


0.26)xl0^ 


2.50 


(1.21± 


0.12)xl0''^ 


(4.57± 


0.46)xlO-' 


3.00 


(8.22± 


0.82)xl0^ 


(7.81± 


0.78)xl0-' 


3.50 


(3.56± 


0.36)xl0^ 


(1.28± 


0.13)xl0^' 


4.00 


(1.15± 


0.12)xl0^ 


(2.04± 


0.20)xl0'^ 


4.50 


(3.03± 


0.30) X 10* 


(3.17± 


0.32) X 10® 


5.00 


(6.89± 


0.69) X 10* 


(4.76± 


0.48) X 10® 


6.00 


(2.60± 


0.26) X 10^ 


(9.52± 


0.95) X 10® 


7.00 


(7.14± 


0.71) X 10^ 


(1.54± 


0.15) X 10^ 


8.00 


(1.50± 


0.15) X 10® 


(2.01± 


0.20) X 10^ 


9.00 


(2.50± 


0.25) X 10® 


(2.18± 


0.22) X 10^ 


10.00 


(3.44± 


0.34) X 10® 


(2.05± 


0.21) X 10^ 
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FIG. 8: The newly derived stellar reactivity of *^Sr(n,p)*^Rb 
is compared to the one given in [s^l (area labelled "new/std"). 
The temperatures Tg are stellar plasma temperatures in GK. 
The shaded area accounts for an error of ±10%. As seen in 
Table 1X1 the experimental data contribute significantly in the 
range 2 < Tg < 4. Also shown is the reactivity from a fit of 
our new result compared to the actual value (curve labelled 
"ncwfit/new"). This shows that the fit accuracy is high. 



—0.562194, ae = 31.1997. The assumed error is con- 
tained in the error given for ag. The coefficients for the 
(p,n) direction are the same, except Cq" = 33.44895;}] ["J g 
and a^*" = -22.3196405. According to Eqs. (gD and ©, 
to obtain the final (p,n) rate the value obtained with the 
seven parameter expression has to be multiplied not only 
by the number densities of the interacting particles but 
also by the ratio of the temperature-dependent partition 
functions of initial and final nucleus 

NpJl'^^{T) = exp 1^1^ + ^ + + a^T^/' + a^T^ 

+ar,Tl'^ + a^XnT^^ , (14) 

^A7^;„(^) = ^^l^Np^n'^^iT) . (15) 

The required partition functions G{T) are provided in 
(lll | as a function of temperature. 

Figure [8] shows a comparison of the new ^^Sr(n,p)^^Rb 
reactivity to the "standard" one of [3^ . At temperatures 
above 3 GK, we see an increase of 10 — 30% compared 
to the previous values. Below 2 GK, the new reactivity 
is 20 — 30% lower than previously. The change in the 
temperature dependence is due to the different proton 
optical potential used. At very low temperature, the re- 
activity becomes less sensitive to the proton potential. 
This explains the ratio becoming almost unity towards 
zero temperature. Also shown in Fig. [5] is a comparison 
between the fit of the new reactivity with the parameters 
above and the reactivity itself. This ratio stays close to 
unity for all temperatures. The deviations between the 
reactivity and its fit are very small and negligible com- 
pared to the other uncertainties involved. 

IV. SUMMARY 

We showed that ~ contrary to common wisdom - a 
large number of endothermic reactions exhibit smaller 
stellar enhancement than their exothermic counterparts 
and are thus preferable for experimental studies. The 
main cause of suppression of the SEF in an endother- 
mic reaction is the Coulomb suppression of transitions 
with low relative energy. This Coulomb suppression of 
the SEF was found to act for reactions with Q < but 
low \Q\ and charged projectiles. Allowing only nucleons, 
a particles, and photons as projectiles or ejectiles, and 
restricting the results to experimentally useful values of 
the SEFs, this effect still appears in 1200 reactions, in- 
cluding a captures relevant in the p proce ss Q and 
proton captures relevant in the rp process [15| and the 
vp process (IQ] . A large number of cases was also found 
for (p,n) reactions which allow to determine astrophysical 
reaction rates relevant to the 7 process @. 

As an example, we measured the astrophysically im- 
portant reaction ^^Rb(p,n)^^Sr close above the thresh- 
old in the energy range relevant for the 7 process. It 
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was shown that in this case it is possible to derive astro- 
physical reaction rates for the (n,p) as well as the (p,n) 
direction directly from the (p,n) data despite of the neg- 
ative reaction Q value. Additionally, our measurement 
confirms a previously derived modification of the global 
proton optical potential used in theoretical predictions. 
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